Charge-density-wave and topological transitions in interacting Haldane model 
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Haldane model [1, is a noninteracting model for spinless fermions showing nontrivial topological 
properties. Effect of the electron-electron interaction on the topological phase pose an intriguing 
question. By means of the Hartree-Fock mean field, the exact diagonalization and the constrained- 
path Monte Carlo methods we mapped out the phase diagram of the interacting Haldane model. 
It is found that interaction breaks down the topological phase and drives the system into the 
charge- density- wave state. Sequence of the two transitions depends on the strength of next-nearest- 
neighbor hopping. Many-body Chern number and the charge excitation gap are used to characterize 
the topological transition. 

PACS numbers: 71.10.Fd, 03.65.Vf, 71.45.Lr ,71.27.+a 



Introduction-Recently, there has been a raising inter- 
est on a class of topological phase of matter: the topo- 
logical insulators, see [2j |3] for reviews. Questions like 
what is the interaction effect on the topological insula- 
tors and how to characterize the topological phase be- 
yond single particle basis pose intriguing challenges, g]- 
[8]. Haldane model [Ij provides a prototype for investigat- 
ing these questions. The model describes non-interacting 
spinless fermions on a honeycomb lattice. What makes 
it interesting is the complex next-nearest-neighbor hop- 
ping ^2, which opens a bulk gap and produces gapless 
edge states. The state is topologically nontrivial and is 
characterized by non-zero Chern number [9j. The topo- 
logical phase breaks down as staggered onsite energy M 
exceeds a critical value. 

The minimal form of interaction in Haldane model is 
the extended Hubbard interaction between the nearest- 
neighbors. The interacting Haldane model reads 

H = Ho^Hi 

Ho = ti^c\cj^it2 ^ijc\cj 

= F^(n,-i)(n,-i) (1) 

The nearest neighbor hopping ti is set as unit in the fol- 
lowing. Uij = ±1 when the hopping path belongs to the 
(anti) clockwise loop. The honeycomb lattice is a bipar- 
tite with two sites (denote as A and B) per unit cell. We 
focus on the half-filling case, i.e. one particle per unit cell 
and the system preserves the particle-hole symmetry. For 
vanishing t2 the extended Hubbard interaction drives the 
system into a charge- density- wave (CDW) state at finite 
value of F [TQHT2 j through a second order phase transi- 
tion. For non- vanishing ^2, there are several questions. 
1. How does t2 affect the CDW transition? What is the 
t2 — V phase diagram looks like ? 



2. Presumably strong enough interaction would breaks 
the topological phase eventually. How does the topolog- 
ical transition organizes itself with the CDW transition? 
Are they separated or merged together? There are sev- 
eral scenarios: i) A continuous or weakly first order CDW 
transition occurs first and there is a finite region where 
small CDW order and the topological phase coexists, ii) 
The two transitions merge: there is a single strongly first 
order transition where the CDW order develops and the 
topological phase breaks at the same time, iii) The topo- 
logical transition occurs first, resulting a topological triv- 
ial phase without CDW long range order. 

3. How to identify the topological transition for an 
interacting system ? Unlike conventional phase tran- 
sitions, the topological trivial/nontrivial phases do not 
correspond to different symmetry groups. The transi- 
tion is of topological nature and can not be described by 
the conventional Ginzburg-Landau- Wilson theory. More- 
over, the model is not exactly solvable once we introduce 
an interaction. 

To address these questions, we perform the Hartree- 
Fock (HF) mean-field, the exact diagonalization 
(ED) and the constrained-path quantum Monte Carlo 
(CPQMC) calculations. A key ingredient in our ED and 
CPQMC calculations is the twisted boundary conditions. 
Average over random twists greatly reduces the finite size 
effect [13HT6]. Moreover, the twisted boundary conditions 
also provide a way of computing many-body Chern num- 
ber [17J: a direct characterization of topological orders. 
Since the topological transition is gap-closing by nature 
[TJ [iHl [19] , the charge excitation gaps are calculated as 
another indicator. 

With these methods we map out the phase diagram of 
CDW and topological transitions and find that HF theory 
predict that a continuous CDW transition occurs before 
the topological phase transition. However many-body 
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FIG. 1: CDW order parameters as a function of V calculated 
by HF (red squares) and GL (green dots) theory for t2 — 
0.5. Blue diamonds denote the charge excitation gap from 
HF calculations. 




FIG. 2: Mean- field t2 — V phase diagrams. Red square line 
defines the CDW phase boundary. Green dot line denotes the 
topological transitions phase boundary, where the CDW order 
parameter exceeds 3v^^2. Inset shows the phase diagram 
taken into account of renormalization of ti (see text). 



calculation reveals that the mean-field calculation misses 
to capture an important possibility that at small t2 the 
topological phase transition happens without developing 
the CDW long range order. 

Hartree-Fock mean- field and Ginzburg- Landau calcula- 
tion-Viider the Hartree-Fock approximation the interac- 
tion term Vfiifij is approximated by {ni)nj + ni{nj) — 
{ni){nj). The many-body Hamiltonian is then written 
into a summation of single particle terms. Introducing 
the CDW order parameter M = zV{{ni^A) — {'riieB))/'^ 
( 2; = 3 is the coordination number), then it couples to 
the density difference of two sub-lattices. The HF mean- 
field Hamiltonian is Fourier transformed into the momen- 
tum space and diagonalized on 200 x 200 meshes of the 
Brillouin zone (BZ). Converged CDW order parameter is 
determined self-consistently. 

For ^2 = we have a second order CDW transition 
at zV = 2.23. Finite next-nearest-neighbor hopping ^2 
pushes the transition point to a larger value of V. It 
is because ^2 term connects sites within the same sub- 
lattices and suppresses the CDW order. However, the 
CDW transition remains to be continuous. Interestingly, 
inside the CDW phase both the order parameter and 
the derivative of the ground state energy shows a kink 
(see Figjl] recognized as the signature of the topological 
transition [20J. Hartree-Fock theory reduces the interac- 
tion problem to the non-interacting Haldane model with 
a staggered onsite energy M. As M exceeds 3V3t2, the 
topological phase transition occurs [Ij. We also calculate 
the charge excitation gap A in the HF energy spectrum. 
It goes to zero right at the topological transition point, 
see Fig{T] The rising of the CDW order parameter and 
the closing of the charge excitation gap divide the phase 
diagram into three regions. The size of the intermedi- 
ate phase enlarges with the increasing of ^2 (see Fig(2|, 



since intuitively, the Haldane phase is more stable with 
a larger value of ^2- 

To further reveal the nature of the transitions, we look 
into the Ginzburg-Landau theory of the model. The par- 
tition function reads Z = J Vc^Vce~^ ^ where the action 

is6'[c'*',c] = drtljl[{dr — IJ^)-\-hk- ^]ipk-^ Hi, a diie Pauli 

matrices. We have introduced ipk = { CkA CkB )^ and 
write Hq on this basis[lj. By performing the Hubbard- 

Stratonovich transformation: e^o "^^^ ^''^"''^^^ = 
/D(^e-^o^^^^^e^>^^^^^^^-^^+^^+^^^ we introduce a 
bosonic field ^ in each unit cell. It acts as the CDW 
order parameter, which couples to the difference of the 
densities on two sub-lattices. Integrating out fermions 
and keeping the zero frequency and momentum compo- 
nent (those contributions dominant over the instability) 
of bosonic field, we have 

s = '-:^(i-zVxM\'+^-^b\<p\\.. (2) 

where X = Efe ^^4rn^ is the CDW susceptibility. Saddle 
point of the effective action of Eqj2] gives the Hartree- 
Fock theory. Second order CDW transition occurs at 
1 — zVx = 0. Near the transition point the order pa- 
rameter follows zVcj) = y\(^zvf^ ' 

In Fig[l|we plot the 
GL prediction of the order parameter. The GL and HF 
result are in good accordance with each other near the 
CDW transition point. However, the kink feature of the 
CDW order parameter is missing in the GL result. It is 
because that GL expansion is only valid near the tran- 
sition point and integration out of the fermions becomes 
invalid at the gap closing point. 

We also like to remark that including the exchange 
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FIG. 3: ED results of Chern number (dots) and charge ex- 
citation gap (triangles) of a 12 sites lattice for different t2. 
20 X 20 twists are used to calculate the Chern number. Inset 
shows the discontinuity of the Chern number point v.s. t2. 

terms V(cJcj)QC^ + h.c. in the mean field decoupling 

amounts to renormalize ti ti + V(cjcj). It renders the 
CDW transition into a first-order transition for t2 < 0.4 
and there is no intermediate phase, see Fig|2] Inset. 

Exact diagonalization and Constrained-path quantum 
Monte Carlo calculations-ln the following, we perform 
ED and CPQMC calculations [21, 22j of the model. In 
CPQMC, the many-body wave function is written as a 
summation of Slater determinants. Imaginary time evo- 
lution operator e~^^ is applied on it to project out the 
ground state. For the present case, auxiliary Ising fields 
are introduced on bonds of the lattice to decouple Hi into 
single-particle form. Summation over these Ising fields 
are performed by Monte Carlo sampling. Different with 
the spinful Kane-Mele-Hubbard model[7l[8], QMC sim- 
ulation encounters sign problems here. Moreover, due 
to the complex next-nearest neighbor hopping t2 and 
the twisted-boundary-conditions, the wavefunctions and 
their overlaps become complex number. The phase prob- 
lem is controlled approximately by the constrained path 
approximation where the overlaps between the Slater de- 
terminants and a given trial wave function are required 
to be positive[23j. Recently, the method has been used 
in the investigation of the spin density wave state [24| [25j 
and ferromagnetism [26j in the Hubbard model. 

Since the twisted boundary condition plays an im- 
portant role in our calculations, we briefly explain it 
here. Under twisted boundary condition, the wave func- 
tion gains a phase when electrons hop around lattice 
boundaries: ^(..., rj + L, ...) = e*^ ®^(..., rj, ...), where 
= {Ox^Oy). The Chern number is calculated by 
C = A. J J d0^d0y[{§^\§i^) - {^|£>], Where |*) is 
the many-body ground state wavefunction with twist 
[17J. Numerical implimentation of above equation fol- 
lows Fukui et al. [18j with one difference: the calculation 




FIG. 4: Comparison of the ED and MC results for a 24 
sites cluster with t2 = 0.1. Inset shows the cluster geome- 
try (PBC). Pink diamonds, green squares and blue triangles 
denote the MC results with A = (noninteracting), A = 1 
(Hartree-Fock) and A = 0.6 reference states. Red line denotes 
ED results on the same cluster. 

was perform on many-body wavefunctions from ED, not 
single particle Bloch states. For other observables like en- 
ergy and CDW structure factor, averaged over random 
twists (O) = / / dOa,dOy{^\0\'^) greatly reduced the fi- 
nite size effect in ED and CPQMC calculations [TSHlS]. 

We first show Chern number versus interaction 
strength V for different t2 in Fig|3] Even for very small 
cluster we see sharp transition of the Chern number 
(from —1 to 0). We see that the topological transition 
point decreases continuously with t2 (Inset of Figj3|. On 
the other hand, it is known that the CDW transition 
point is finite even for vanishing t2 [TT| . There- 
fore for small enough t2 the topological transition would 
occur before the CDW transition. This result is not 
captured in the mean-field phase diagram Fig|2] since 
there the only way of breaking the topological phase 
is to develop the CDW long range order. However in 
fact for small t2 the Haldane phase is so fragile that 
short range CDW fluctuations may destroy it, resulting 
a gapped phase with neither topological nor CDW long 
range order. We also calculate the charge excitation gaps 
A = {En^i + En_i — 2En)/2, where En is the lowest 
energy in the half-filling sector, En±i is the lowest en- 
ergy in the sector with one more (less) particle. Since the 
particle-hole symmetry is respected, the above equation 
can be simplified to A = {En-i + zV/2 — En). From 
Figj3]we see that the minimum point of the PBC gap is 
almost coincide with the jumping point of Chern num- 
ber. Thus, we link the PBC gap-closing point with the 
topological transition. Since Chern number calculations 
with ED do not scale up to large clusters, in the following 
we calculate the charge excitation gap with CPQMC for 
the signature of topological transition. 

Figj4] shows the benchmarks of the CPQMC results 
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FIG. 5: TABC CDW structure factor (dots) and PBC charge 
excitation gaps (triangles) vs. interaction strength V for t2 = 
0.1 (upper panel) and t2 = 1.0 (lower panel). Red and green 
line indicates the CDW and topological phase transition point 
from finite size scaling. Insets show a zoom in view of the 
charge excitation gap for t2 = 0.1 and cluster geometries. 



FIG. 6: Phase diagram determined by CPQMC calculations. 
Red squares and green dots denote the phase boundaries of 
the CDW and topological transitions. Error bars indicates the 
uncertainties of determining the phase transition boundaries. 
Lines are guides to the eye. 



with the ED results. If PBC is adopted and the number 
of columns is a multiply of 3 then K-point is a validate 
point in the momentum space. The reference state for 
CPQMC calculation is obtained from the HE calcula- 
tion with V = XV where A is an interpolation parame- 
ter, resembles the tuning parameter in hybrid function- 
als. The noninter acting state corresponds to A = and 
the "true HE state" corresponds to A = 1. We sample 
the ground state energy and the CDW structure factor 
^CDW ^ i_ ^.^^.((nf - nf){nf - nf)) in MC calcula- 
tions. The two MC energies crosse at V = 3.0. On both 
sides, we take the physical quantities from the one with 
lower energies. Practically, we also adopt the A = 0.8 for 
t2 = 1.0 and A = 0.6 for t2 = 0.1. Both approaches give 
results show good accordance with exact ones, see Figj4] 
For the 24 sites lattice, the CDW structure factor shows 
a jump dit V = 3.0. 

We show the CDW structure factor and the charge ex- 
citation gap from CPQMC calculations in Fig|5] The 
CDW structure factors are continuous and obey the scal- 
ing law of second order transition = /[(^ ~ K)^'^]? 
with / a universal function independent of system size. 
By fitting the scaling law we extract the CDW transi- 
tion point Vc. Gap closing is interpreted as signature of 
the topological transition. We mark the CDW and topo- 
logical transition point in the figure. Widths of them 
indicates uncertainty of the fitting. Even with the uncer- 
tainties, we could clearly see the topological transition 
occurs before the CDW transition for ^2 = 0.1 and after 
it for t2 = 1.0. 

We conclude our many-body calculations with phase 
diagram shown in FigjOj Different with the mean-field 
prediction in Figj2j the CDW and topological transition 



lines intersect. Above the intersection, the nature of two 
transitions is mean-field like, i.e. long-range CDW order 
develops and then destroys the Haldane phase. Albeit the 
two phase boundaries shift against the mean-field values. 
Below the intersection the topological phase is destroyed 
before the CDW order develops. As ^2 ^ the topo- 
logical boundary goes to zero while the CDW transition 
point remains at a finite value. 

Outlook-First^ to better characterize the topological 
transition in an interacting system, we will look into the 
calculation of many-body Chern number or entanglement 
entropy /spectrum within CPQMC method. Second, the 
nature of topological transition in small t2 region de- 
serves further study. Is the topological phase been de- 
stroyed purely due to fluctuation effect or instability to- 
ward other order, say Kekule distortion? Third, results 
reported here may shed light on phase diagram of the 
Kane-Mele-Hubbard model where Z2 topological order 
competes with antiferromagnetic order. 
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